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Exclusive nonleptonic bottom meson decays are studied in the covariant osillator quark 
model using the factorization assumption. The main feature of this model is that it can 
simultaniously be applied to both heavy heavy and heavy light transitions, satisfying 
the constraints of the heavy quark effective theory (HQET) in the appropriate limit. The 
results obtained are in overall agreement with the present experimental data for various B 
decays. 
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^ ' SI. Introduction 
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The description of exclusive nonleptonic decays of B mesons represents an im- 
portant and complicated theoretical problem. These decays are nonperturbative in 
nature and cannot be calculated reliably from the QCD Lagrangian. In contrast 
, to exclusive semileptonic decays, where the weak current matrix elements between 

meson states are involved, nonleptonic decays require the evaluation of hadronic 
matrix elements of four local quark operators. To simplify the analysis it is usually 
assumed that the matrix element of the current-current weak interaction factorizes 
I into the product of two single current matrix elements. Thus the problem reduces 

to the calculation of the meson form factors, which are contained in the hadronic 
matrix elements of weak currents as in the case of semileptonic decays, and of the 
meson decay constants, describing the leptonic decays. This makes the factorization 
hypothesis a very appealing assumption. Although it is very difficult to prove 
the factorization hypothesis theoretically within our present understanding of QCD 
nonperturbative effects, this hypothesis is expected to be valid to a rather good 
approximation in the case of transitions with large energy release, such as heavy 
B decays, since the final mesons carrying large momenta escape from the region of 
interaction, thereby minimizing the effects of a final state interaction. Several tests 
have been made to prove its validity phenomenologically, and, it has been shown to 
work well for the description of B meson decays into a D or D* and a light meson. 
Once the factorization assumption is made, nonleptonic decays are related to the 
corresponding semileptonic decays. 
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In this paper we wish to calculate the branching ratios of the exclusive nonlep- 
tonic two-meson decays of B^, B~ and Bs mesons in the framework of the covariant 
oscillator quark model (COQM) on the basis of factorization. One of the most im- 
portant motivations of the covariant oscillator quark model (COQM) is to describe 
covariantly the centre of mass motion of hadrons, retaining the considerable sucesses 
of the non-relativistic quark model with regard to the static properties of hadrons. 
A key element of the COQM for achieving this is its direct treatment of the squared 
masses of hadrons, in contrast to the mass itself as done in conventional approaches. 
This makes the covariant treatment simple. The COQM has been applied to various 
problems^) with satisfactory results. Recently, Ishida et al have studied the semi- 
leptonic weak decays of heavy hadrons using this model and derived the same 
relations of weak form factors as in HQET.*^) Furthermore, the predicted spectra 
for B {D*,D)lv were shown to fit experimental data quite well. Keeping this 
success in mind, we extend application of the COQM to the nonleptonic decays of 
5°, B- and Bs mesons. 

This paper is organized as follows. In §2 we present the expressions for nonlep- 
tonic decay amplitudes in the factorization approximation. In §3 we present a brief 
description of the covariant oscillator quark model. Using this model we have eval- 
uated the form factors and obtained the decay rates for various nonleptonic decay 
processes. The decay modes B —> D*V are considered in §4. Section 5 contains our 
results and discussion. 

§2. General formalism 

Neglecting the penguin contribution, the effective Hamiltonian describing the 
decays under consideration is given by 

T^eff = ^ Kb Vq,q. [Ci(m6)Oi + C2(m6)02] (2-1) 

with 

Oi = {mjY{cb)n and O2 = mY {Mj) n . (2-2) 

where Ci and C2 are the Wilson coefficients, and the quark current {qiqj)fi denotes 
the usual {V—A) current, qi and qj are two types of quark flavors that are hadronized 
to the P or y mesons. 

The factorization approach to two-body nonleptonic decays B DM implies 
that the decay amplitude can be expressed by the product of one particle matrix 
elements: 

a^{D\{ch)^\B){M\{qiq^r\Q) 



{DM\H,s\B) = -kVcbV^ii, 



V2 



+a2{M\{qib)^\B){D\{cqjr\0) 



(2-3) 
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Here ai = Ci + C^j^c and 02 = C2 + CxjNci where Nf. represents the number of 
colors. 

It should be noted that in writing Eq. (2-3) wc have discarded the contribution 
of color octet currents which emerges after the Fierz rearrangement of color singlet 
operators. Clearly these currents violate factorization since they cannot allow tran- 
sition to the vacuum states. In the factorization approximation one can distinguish*^ 
three classes of B meson decays : the 'class I' transitions such as B^ M2 , 

where only the term ai contributes (both mesons produced by charged currents ); 
'class ir transitions, such as B^ — > Mj^M^j, where only the term 02 contributes 
(both mesons are produced by neutral currents); and 'class III' transitions, such as 
B~ — ^ M1M2 , where both terms contribute coherently. 

In order to evaluate the transition amplitudes we use the following matrix ele- 
ments: 

(P(p)|(g,g,nO) = -i/pp^ 
(nP,e)lfer|0) = My /ye^ 

{ai{p,e)\{mjY\Q) = Ma, fa, e^ (2-4) 

where P, V and ai represent the pscudoscalar, vector and the axial vector mesons, 
respectively. To evaluate the hadronic form factors we use the COQM. These are 
presented in the next section. 

§3. Model Framework, hadronic form factors 
and decay width of PP, PV and B^ VP 

The general treatment of COQM may be called the "boosted LS-coupling scheme," 
and the wavefunctions being tensors in ?7(4) x 0(3, l)-space, reduce to those in 
5'C/(2)spin X 0(3)orbit-space in the nonrelativistic quark model in the hadron rest 
frame. The spinor and space-time portion of the wave functions separately satisfy 
the respective covariant equations, the Bargmann-Wigncr (BW) equation for the 
former and the covariant oscillator equation for the latter. The form of the meson 
wave function has been determined completely through the analysis of mass spectra. 

In COQM, the meson states are described by bi-local fields <^^(xi^, X2^), where 
xiij,{x2p) is the space time coordinate of the constituent quark (antiquark), A = 
(a, a) (B = (b, /?)) describing its flavor and covariant spinor. Here we write only the 
positive frequency part of the relevant ground state fields: 

^2/.) = e^''-^ U{P)i fabix^; P) , (3-1) 

where U and / are the covariant spinor and internal space-time wave functions 
respectively, satisfying the Bargmann-Wigner and oscillator wave equations. The 
quantity x^{X^) is the relative (CM) coordinate, = xi^ — X2ij,{Xn = mixi^ + 
'm2X2fj) / {mi + m2), where the mi represent the quark masses). The function U is 

*^ The contributions due to quark annihilation processes, which are expected to be small, are 
neglected in the calculation of this paper. 
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given by 

UiP) = ^ [{-l5Ps{v) + ilM^m + ■ 7)] , (3-2) 

where Ps{Vs) represents the pseudoscalar (vector) meson field, and = P^/M 
[P^,(M) is the four momentum (mass) of the meson]. The function U, being repre- 
sented by the direct product of quark and antiquark Dirac spinors with the meson 
velocity, is reduced to the non-relativistic Pauli-spin function in the meson rest frame. 
The function / is given*) by 

/(x,;P) = ^exp(-f (4 + 2!^)) . (3-3) 

The value of the parameter P is determined from the mass spectra^-* as /^Tr/p/ai = 
0.14, Pk/k* =0.142, /3^/^.=0.148, Pd, = 0.154, Pb =0.151 and =0.160 (in units 
of GeV^ ). 

The effective action for weak interactions of mesons with VF-bosons is given by 

Sw = j d*xid'^X2{^F,P'{xi,X2)i'y^{l + ^^)^i,p{xi,X2))W^^q{xi) , (3-4) 

where we have denoted the interacting (spectator) quarks as 1 (2). The CKM matrix 
elements and the coupling constant are omitted. This equation is obtained from 
consideration of Lorentz covariancc, assuming a quark CTirrcnt with the standard 
V — A form. In Eq. (3-4), '^i^p {^f,p') denotes the initial (final) meson with definite 
four momentum P^^P'^), and is the momentum of W boson. The function 3 is 
defined by ^ = —74(^^74, and ( ) represents the trace of Dirac spinor indices. Our 
relevant effective current Jn{X)pi^p is obtained by identifying the above equation 
with 

Sw = l d^XJ^{X)p,,p W^{X)g . (3-5) 
Then Jij,{X = 0)p'^p = J^, is explicitly given as^^ 

= I'i\w)^MM' X [Ps{v')Ps{v){v + v')^ 

+ y\{v')Ps{v){e^xai3v'aVi3 " h^i{w + 1) - v^v'^] , (3-6) 

where M{M') denotes the physical masses of the initial (final) mesons. It should be 
noted that in the pure confining limit, the masses of the ground state mesons are 
equal to the simple sums of their constituents, which are much different from the 
physical masses in the case of light quark pseudoscalar mesons, such as tt and K. 
Therefore we do not consider the transitions -B — >■ tt and B ^ K in our analysis as 
the rcliabilty of the results is less for these transitions. The quantity I'^^{w), which is 

*^ In this paper we employ the pure-confining approximation, neglecting the eflFect of the one- 
gluon-exchange potential, which is expected to be good for the heavy /light-quark meson system. 
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the overlapping of the initial and final wave functions, represents the universal form 
factor. It describes the confined effects of quarks and is given by 

= ^fn, 2^ eM-G{w)) ; C{w) = (p - p'f + App'w' , (3-7) 



and 



G{w) = 



mt 



2C{w) 



.r, I / M \2 /M'\^ MM' 



\M'J 



(3-8) 



where Ms(M^) represents the sum of the constituent quark masses of the initial 
(final) meson, and m„ is the spectator quark mass. 

The form factor function I'*{w) for B D{D*)decays corresponds to the Isgur- 
Wisc function ^{w) in HQET. 6) At the zero recoil point w = 1, the value of I'^^{w) 
is given by 

''V-D-pf^. (3-9) 

In the heavy quark symmetry limit (3 = P' , so Eq. (3-9) correctly reproduces'''^ 
the normalization condition of HQET, i.e., ^(w = 1) = 1. However, HQET, as it 
is, predicts nothing about the Isgur-Wisc function except for the zero recoil point, 
while in COQM the form factor functions can be derived at any kinematical point of 
interest due to the fact that the center of mass motion of the meson there is treated 
covariantly, as was mentioned in §1. In addition, the COQM form factor I'^^{w) 
is also applicable for the heavy-to-light transition processes, while HQET does not 
provide anything for this sector. 

After obtaining the effective current in the COQM, the decay widths for various 
B ^ D and B ^ D* decay modes can be obtained with Eqs. (2-3), (2-4) and (3-6). 
These are as follows: 



nB{v) ^ D{vi)P{v2)) = iVcbV,,,/ |p 



B 



ai fp VMBMD{lf{wi)){l + wi) {Mb - Md) 



1 2 



+ as fo VMBMp{lf{w2)){l + W2) {Mb - Mp) 



(3-10) 



r{B{v)^D*{vi)P{v2)) 



167rM| 
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ai ff 



+ a2 fD* J^(^fM) (Mb + Mp) 



(3-11) 



r{B{v) ^ D{vi)V{v2)) 
X ai fv 



^2 

F W T/ |2 |„|3 



+ a2 fD J^ilfM) (Mb + Mv) 



(3-12) 



where we have taken wi = v ■ vi and W2 = v ■ V2- Here |p| is the cm. momentum of 
the emitted particles in the rest frame of the B meson. As stated earher, only the 
ai term contributes to class I decays, and only the 02 term contributes to class II 
decays, while both ai and 02 terms contribute coherently to class III decays. 

The COQM is also applicable to heavy to light transitions, as well as heavy 
to heavy transitions, such as B ^ p and B K*. In this case the above for- 
mulas are changed by replacing {D*{vi), Iq''{wi),Vcb) in B ^ D* transition with 
{p{vi),lf{wi),Vy,b) and {K*{vi),lf{wi),Vub) in B ^ p and B ^ K* transitions, 
respectively. 

§4. Decay rate, polarization and angular correlation in the decays 

B^VV 



The helicity amplitude for the decay process B{p) — > D*{ki,ei)V{k2,e2) can be 
expressed by three invariant amplitudes, a, b and c. It is given following Ref. 10) as 



Hx = e;„' e 



(A)* (A)* 
In *^2v 



ag^'' + 



b I, ,j ic 



Md*Mv 



Md*Mv 



laPl3 



(4-1) 



The coefficients a, b and c describe the S-, P- and D- wave contributions to the two 
final vector particles. In the present framework these are given as 

a = ^V^Vq,q.[alfvMv^/M^M^* I'*{wi)(l + w^) 



Of 



+ a2fD*MD*VMBMv I'^''{W2){1+W2)], 

aifvM^ 



aifvMir 



2 r\^,) + a2fD*Ml, 



Mb 



I My 

Mb 
Mb 



W2 



W2] 



■ (4-2) 
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The helicity amplitudes are given as 
H±i = a lb ysP^—l c 



where x is defined by 



and obeys 



X = 



and Hq = —ax — b{x^ — 1) 

ki-k2 _ M| - M|,. - M2 



Md*Mv 



x^ = 1 + 



2Md*Mv 



The corresponding decay rate can be obtained as 



r{B D*V) = 



2|ap + \xa + {x"^ - 1)6^ + 2{x^ - l)\c\ 



(4-3) 

(4-4) 
(4-5) 

(4-6) 



The decay distribution is parametrized by the coefficients 



r 

oil = 
a2 = 



\H. 



+1 



I' + l^-iP 



|i7oP + |^^+iP + |i/-iP' 

Re {H+iH^ + H^iH^) 
|i7oP + |i?+iP + |i?-iP' 

Re(ff+i^-i) 
\Ho\^ + \H+i\^ + \H^i\^' 



r 



\Ho\'' + \H+i\^ + \H^i\ 

Im (iJ+iiJo* - g-iffp*) 
|i7oP + |i^+iP + |i/-iP 



P2 



\Ho\'' + \H+i\^ + \H^i\ 



(4-7) 



In general, the dominant terms in the angular correlations are Ft/F, Fl/F, a\ 
and a2 ■ The terms (Si and (32 are small since they are nonvanishing only if the helicity 
amplitudes -ff+i, -ff-i and Hq or the invariant amplitudes a, h and c, respectively, 
have different phases. 

In the case of heavy to light transitions, B ^ p and B K* , the corresponding 
formulas are obtained by the procedure explained as the end of the last section. 



§5. Results and conclusion 



In order to make the numerical estimate, we use the following values of various 
quantities. The quark masses (in GeV) are taken as niu = ma = 0.4, = 0.55, 
rric = 1.7 and = 5.^^ The particle masses and lifetimes are taken from Ref. 15). 
The relevant CKM parameter values used are Vcb = 0.0395, Vcs = 1-04, Vcd = 0.224, 
Vy,d = 0.974, Vus = 0.2196 and Vub = 0.08 x Vcb = 0.00316. The decay constants are 
taken as U = 130.7, fx = 159.8; fx* = 214; fp = 210; ^ 220, fo* = 230, 

/d, = 240, fo* = 260^3) and fa, = 205 1'') (in MeV). The decay constant fj/^ is 
determined from the value of F^J/tp e'^e~): 

fj/^ = Y i {-^) ^(-^/^ ^ e+e-)Mj/^ = 404.5 MeV . (5-1) 
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Tabic I. Branching ratios of nonleptonic B decays in the COQM. Note that the values do not 
include the possible contribution from the penguin diagram, which is generally expected to 
be of order 10~® ~ 10~^. It is shown that the penguin diagram does not contribute to 
p+ , p+ Dt' , K*° D° and K*°D*°. The contribution to p+n'^p+p- and p+a^ is about 
3 X 10"''.^^ 



modes 


This work 


Expt. [15] 


Class I 




3.00 xlO" 


■3 


(3.0 ±0.4) X 10"^ 


D+K- 


2.28 X 10" 


-4 


- 


D+D- 


4.345 X 10 


-4 


- 


d+d: 


10.9 X 10" 


-3 


(8.0 ±3.0) X 10"^ 


D+p- 


7.406 X 10 


-3 


(7.9 ± 1.4) X 10"^ 


D+K*- 


3.84 X 10" 


-4 


- 


D+D*- 


3.29 xlO" 


■4 


< 1.2 X 10"^ 


D+or 


0.86 x 10" 


-2 


(1.0 ±0.5) X 10"^ 




6.51 x 10" 


-3 


(6.0 ±3.3) X 10"^ 


D'+n- 


3.07 X 10" 


-3 


(2.76 ± 0.21) X 10"^ 


D*+K- 


2.28 X 10" 


-4 




D*+D- 


3.14 X 10" 


-4 


< 1.8 X 10"=* 


D*+DJ 


7.615 X 10 


-3 


(9.6 ±3.4) X 10"^ 


D*+p- 


8.91 X 10" 


-3 


(6.7 ±3.3) X 10"^ 


D*+K'- 


4.89 X 10" 


-4 




D*+D*- 


8.74 X 10" 


-4 


< 2.2 X 10"^ 


D*+Dt- 


2.49 X 10" 


-2 


(2.0 ±0.7) X 10-2 


D*+a:[ 


0.99 X 10" 


-2 


(1.30 ± 0.27) X 10-2 




2.17 X 10" 


-5 


< 7 X 10-" 




4.63 X 10" 


-5 


< 8 X 10"'' 


P^TK 


6.53 X 10" 


-6 


< 8.8 X lO-'^ 


P+P~ 


1.83 X 10" 


-5 


< 2.2 X 10-^ 


p+a~ 


1.94 X 10" 


-5 


< 3.4 X 10-^ 


Class II 




0.649 X 10 


-4 


< 3.9 X 10-* 


D*«p" 


1.22 X 10" 


-4 


< 5.6 X 10-'' 




1.504 X 10 


-3 


(1.35 ±0.18) X 10-^ 




1.16 X 10" 


-6 




j~^*0 


2.25 X 10" 


-6 





The parameters ai and 02 appearing in these decays, which have recently been 
determined from the CLEO data, are ai = 1.02 and 02 = 0.23. i'^) Using these 
values we obtain the branching ratios for B^, B~ and Bg mesons, which are tabulated 
in Tables I, II and III, respectively. The overall agreement between the predicted and 
experimental data is quite remarkable. Here, it may be worthwhile to note that the 
relevant processes are relativistic and the form factor functions /(cu) play significant 
roles: For example, the velocity of the final D'^ in the process B^ — > D^tt^{B^ — ^ 
D+D*-) is VD = 0.78c(0.70c), and the value of / is / = 0.53(0.65). The polarization 
and the angular distribution parameters for B D*V decay modes arc presented 
in Table IV. The polarization fraction (ri/r) for B° D*+ p' (0.883) agrees weU 
with experimental value 0.93 ± 0.05 ± 0.05. For the decay mode B° K*^JI'^, 
FL/r (0.428) also agrees with the recent CLEO data 0.52 ± 0.07 ± 0.04. 



Non-Leptonic Two-Meson Decays of B-Mesons 



9 



Tabic II. Branching ratios of nonleptonic decays in the COQM. The numbers for _D"aj^, 
D*" ,D'^ ■K~ ,D*^ ■K~ ,0^^ K~ ,D*^ K~ , p°7r~ and p^a^ in class III do not contain the contri- 
butions from the color-suppressed a2-term in Eq. (2.3). Note also that the values do not 
include the possible contribution from the penguin diagram, which are generally expected 
to be of order 10~^ ^ 10"^. It is shown that the penguin diagram docs not contribute to 
p'^ D~ , p'^ Dl~ , p~ J / ip , K*~ and K*~ D*". The contribution to p^v~ ,p^p~ and p°a]~ is about 
3 X 10"'^.^^ 



B modes 


This work 


Expt. [15] 


Class I 


D^D- 


4.53 xlO" 


■3 




D°D- 


1.16 X 10" 


-2 


(1.3 ±0.4) X 10"^ 


D°D*- 


3.48 X 10" 


-4 




D°D*- 


9.07 X 10" 


-3 


(9.0 ± 4.0) X 10-3 


D*°D- 


3.34 X 10" 


-4 




d*'>d: 


0.81 xlO" 


-2 


(1.2 ±0.5) X 10"^ 


D'^D*- 


9.23 X 10" 


-4 






2.625 X 10 


-2 


(2.7 ±1.0) X 10"^ 


p'DT 


1.15 xlO" 


■5 


< 4 X 10"* 


p'DT 


2.45 X 10" 


-6 


< 5 X 10"* 


Class II 


K'-J/i, 


1.587 X 10- 


-3 


(1.47 ± 0.27) X 10"^ 


p-Jjib 


5.79 X 10" 


-6 


< 7.7 X 10"* 


K*-& 


1.23 X 10" 


-6 




K*-D*° 


2.38 X 10" 


-6 




Class III 


^v- 


1.004 X 10 


-2 


(1.34 ±0.18) X 10"'' 


D'"p- 


1.264 X 10' 


-2 


(1.55 ±0.31) X 10"^ 


D°a^ 


6.89 X 10" 


-3 


(5.0 ±4.0) X 10"3 




1.04 X 10" 


-2 


(1.9 ±0.5) X 10"^ 


D\- 


3.18 X 10" 


-3 


(5.3 ±0.5) X 10"^ 




3.25 X 10" 


-3 


(4.6 ± 0.4) X 10"=^ 


D°K- 


2.41 X 10" 


-4 




D°K*- 


5.34 X 10" 


-4 






2.42 X 10" 


-4 




D*°K*- 


7.13 X 10" 


-4 




p%- 


3.45 X 10" 


-6 


< 4.3 X 10"^ 


p°ar 


1.03 X 10" 


-6 


< 6.2 X 10"* 


p'p- 


1.45 X 10" 


-6 


< 1.0 X 10"^ 



In this paper we have calculated the branching ratios of the exclusive nonlep- 
tonic two-meson decay of B mesons using the covariant oscillator quark model based 
on the factorization approximation. The applied form factors are consistent with 
the predictions of heavy quark symmetry. The overall agreement of our predictions 
for two meson nonleptonic decays of B mesons with the existing experimental data 
suggests that the factorization approximation works well and the estimation of con- 
finement effects by the COQM is valid. 
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Table III. Branching ratios of nonleptonic decays in the COQM model. 



modes This work Expt. [15] 





0.248 xlO"^ 


< 13% 


D+K- 


1.89 X 10"^ 


_ 


D+D- 


3.73 X 10"* 


_ 


D+D7 


9.57 X 10"^ 


_ 


D+p- 


6.15 X 10"^ 




DfK*- 


3.196 X 10"'' 




D+D*- 


2.86 xlO"" 






7.48 x 10"^ 




D+a^ 


4.86 X 10"'' 




D*+n- 


2.55 X 10"^ 




D*+K- 


1.9 X 10"* 




D*+D- 


2.735 X 10"'' 




Dl+D- 


6.68 X 10"'^ 




Dl+p- 


7.49 X 10"^ 




D*+K*- 


4.1 X 10"" 




d:+d*- 


6.85 X 10"* 





Dt+D*- 2.216 X 10"^ 
Df^a^ 8.39 X 10"^ 



Table IV. Polarization and angular correlation parameter oi B ^ D*V decays. 



Decay modes 


rL/r 


Ql 


Q2 


^ D'+p- 


0.883 


-0.416 


0.04 


B° -c D'+D*- 


0.538 


-0.662 


0.176 


B° D*+D*s- 


0.515 


-0.665 


0.187 




0.428 


-0.605 


0.142 


B- -> D*^^p- 


0.855 


-0.446 


0.044 


B- D'°D'- 


0.538 


-0.661 


0.176 


B- D*°D*- 


0.515 


-0.665 


0.187 


B- ^K*-J/i) 


0.428 


-0.605 


0.141 
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